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We provide a general dynamical approach for the quantum Zeno and anti-Zeno effects in an 
open quantum system under repeated non-demolition measurements. In our approach the repeated 
measurements are described by a general dynamical model without the wave function collapse pos- 
tulation. Based on that model, we further study both the short-time and long-time evolutions of the 
open quantum system under repeated non-demolition measurements, and derive the measurement- 
modified decay rates of the excited state. In the cases with frequent ideal measurements at zero- 
temperature, we re-obtain the same decay rate as that from the wave function collapse postulation 
(Nature 405, 546 (2000)). The correction to the ideal decay rate is also obtained under the non- 
ideal measurements. Especially, we find that the quantum Zeno and anti-Zeno effects are possibly 
enhanced by the non-ideal natures of measurements. For the open system under measurements with 
arbitrary period, we generally derive the rate equation for the long-time evolution for the cases with 
arbitrary temperature and noise spectrum, and show that in the long-time evolution the noise spec- 
trum is effectively tuned by the repeated measurements. Our approach is also able to describe the 
quantum Zeno and anti-Zeno effects given by the phase modulation pulses, as well as the relevant 
quantum control schemes. 
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I. INTRODUCTION 

Quantum Zeno effect (QZE) [l[ and quantum anti- 
Zeno effect (QAZE) are among the most interesting 
results given by quantum mechanics. The two effects de- 
scribe the evolution of a quantum system under repeated 
(or continuous) measurements. QZE shows that, for a 
closed quantum system with finite-dimensional Hilbcrt 
space, the unitary evolution can be inhibited by the re- 
peated measurements. Not surprisingly, the similar inhi- 
bition can also occurs in the Rabi oscillation of a quantum 
open system coupled with an environment. 

On the other hand, for an unstable state in an open 
quantum system with dissipation, QZE shows that the 
time-irreversible decay of such a state can be suppressed 
when the quantum measurements are frequent enough. 
In addition to the QZE in such an extreme limit, the de- 
cay rate of the unstable state in a dissipative system is 
also possible to be increased if the measurements are re- 
peated with an "intermediate" frequency. That is known 
as quantum anti-Zeno effect. QZE and QAZE have at- 
tracted much attention since they were proposed. So 
far the two effects have been experimentally observed in 
the systems of trapped ions [1, [f| , ultracold atoms 0- 
[U, molecules H and cavity quantum electric dynamics 
(cavity QED) jlo|. 

The QZE and QAZE were previously derived with the 
postulation of wave packet collapse in the quantum mea- 
surement. After the initial proposals of the two effects, 
many authors have discussed the dynamical explanation 



of QZE and QAZE via the unitary quantum mechanical 
approach without wave packet collapse. For the unitary 
evolution of a closed system or the Rabi oscillation of an 
open system, many authors pointed out that the QZE 
in these s yste ms can also be explained in a dynamical 
approach |lll - l33 |. These explanations are based on the 
considerations of the quantum dynamics of the quantum 
system as well as the apparatuses during the measure- 
ments. In these complete dynamical analysis, QZE can 
appear naturally during the evolution of the reduced den- 
sity matrix of the system under the frequent interactions 
with the apparatuses, and the postulation of wave packet 
collapse is not required. These explanations of QZE were 
previously given for many special systems and set-ups 
for measurements (e.g., two- level atoms under the detec- 
tion with laser pulses). Recently some of us (D. Z. Xu, 
Qing. Ai and C. P. Sun) provided a general dynami- 
cal proof of QZE for any closed system under quantum 
non-demolition measurements 
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For the QZE and QAZE of the decay of an unstable 
state in a quantum open system with dissipation, there 
are also similar discussions in which the d yna mical pro- 
cess of the measurements is included |32|-;38j . In Refs. 
[3~4T - |37| , where the quantum open system is also assumed 
to be under repeated measurements, the dynamics of ev- 
ery individual measurement is discussed in detail. Never- 
theless, in these references the total survival probability 
of the initial state of the open system is just intuitively 
calculated as the product of the one after each measure- 
ment, rather than derived from the dynamical equation 
for the total evolution process. This treatment is appro- 
priate for the cases with ideal projective measurements. 
However, as we will show in Sec. Ill, in the cases of 
general measurements which could be non-ideal, a more 
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first-principle analysis is required. There have also been 
some full dynamical discussions of QZE and QAZE in 
the limit of continuous measurement [HI, [33] (i.e., the 
cases where the system keeps interacting with the mea- 
surement apparatus during the total evolution time, or 
the time interval between two measurements are much 
smaller than the duration time of each measurement) 
or the case of repeated measurements [38| . However, in 
these discussions the measurements are performed via 
several special physical systems, rather than general ap- 
paratus. To our knowledge, so far there has not been a 
general dynamical explaination for QZE and QAZE of 
the unstable states in a quantum open system for arbi- 
trary environment spectrum, environment temperature 
and system-apparatus coupling. 

In this paper, we provide a general dynamical approach 
for the QZE and QAZE of unstable states in a quantum 
open system. For simplicity, we illustrate our central 
ideas with a two-level system (TLS) coupled to a heat 
bath of a multi-mode bosonic field. As pointed in Sec. 
VI, it is straightforward to apply our model in open sys- 
tems and environments of other kinds. 

In our approach, every measurement process is de- 
scribed by the dynamical coupling between the quantum 
open system and the apparatus. The postulation of wave 
function collapse is not used. To describe the repeated 
measurements, we use a multi-apparatus model that in 
each measurement, the open system interacts with an in- 
dividual apparatus [HJ , while all the other apparatus are 
left alone. As proved in Sec. II, regarding to the evolu- 
tion of the density operator of the quantum open system, 
this multi-apparatus model of repeated measurements is 
equivalent with the intuitive model where all the mea- 
surements are given by the same apparatus and the state 
of this apparatus is initialized before each measurement. 
We also point out, in our model there is no special re- 
striction on the details of the measurement process, e.g., 
the coupling between the to-be-measured system and the 
heat bath. To be more practical, we concentrate on the 
case of " repeated measurements," i.e., the duration time 
of each measurement is much smaller than the time in- 
terval between two measurements. 

Using the multi-apparatus dynamical model of re- 
peated measurements, we completely investigate the ef- 
fect of the measurements on both the short-time and 
long-time evolution of the TLS, and derive the QZE and 
QAZE in various cases. For the short-time evolution, 
with time-dependent perturbation theory, we obtain the 
short-time decay rate of the excited state of the TLS, and 
show that when the TLS is under repeated ideal projec- 
tive measurements, the short-time decay rate given by 
our dynamical model is exactly same as the one from 
the wave packet collapse postulation Q. In the case of 
imperfect measurements, the corrections due to the non- 
ideality and the finite duration time of the measurements 
can also be naturally obtained. We prove that in the case 
of non-ideal measurements, QZE also occurs when the 
measurements arc frequent enough. 



For the long-time evolution of the TLS, we derive the 
rate equation of the TLS by calculating the complete 
long-time evolution of the TLS and the environment, 
rather than by resetting the state of the environment 
after each measurement. Since the Markovian approxi- 
mation is not applied, our rate equation can be used for 
the cases of heat baths with arbitrary correlation time 
and the cases of measurements repeated with arbitrary 
frequency. The rate equation clearly shows that the spec- 
trum of the noise is effectively tuned by the repeated mea- 
surements. Especially, we show that when the measure- 
ments are frequent enough, the time-local approximation 
and coarse-grained approximation will be applicable. In 
this case, at zero temperature, the decay rate given by the 
long-time rate equation is the same as the short-time de- 
cay rate from our short-time perturbative calculation in 
Sec. IV. Furthermore, as a result of the counter-rotating 
terms in the coupling Hamiltonian of the TLS and the 
heat both, the decay rate of the ground state of the TLS 
may be varied to a non-zero value by the periodic mea- 
surements, even at zero temperature. 

Our model is also able to describe the QZE or QAZE 
given by periodic phase modulation pulses [39144 1| rather 
than measurements, since the former one can be con- 
sidered as a special kind of non-ideal measurement with 
complex decoherence factor. Our calculation also pro- 
vides a microscopic or full-quantum theory for the re- 
cent prop osals of stochastic control of quantum coher- 
ence j42H44l |. where the repeated projective measure- 
ments are semi-classically described as a stochastic term 
in the Hamiltonian. 

This paper is organized as follows. In Sec. II we post 
our multi-apparatus dynamical model for the repeated 
quantum non-demolition measurements. In Sec. Ill we 
show our dynamical description for the dissipative TLS 
under repeated measurements. In Sec. IV we calculate 
the short-time decay rate of the excited state of the TLS 
via time-dependent perturbation theory, and discuss the 
effects given by non-ideal measurements or phase modu- 
lation pulses to QZE and QAZE. In Sec. V we consider 
the long-time evolution of the system, derive the effective 
noise spectrum experienced by the TLS and obtain the 
rate equation for the long-time evolution at finite tem- 
perature. There are some conclusions and discussions in 
Sec. VI. 



II. MULTI-APPARATUSES FOR REPEATED 
MEASUREMENTS 

In this paper, we discuss the quantum evolution of a 
dissipative TLS under repeated quantum non-demolition 
measurements. There are two possible models to describe 
the repeated quantum measurements, i.e., the single- 
apparatus model and the multi-apparatus model. In the 
single-apparatus model, all the measurements arc com- 
pleted via the same apparatus whose state is "initial- 
ized" to a special one in the beginning of each measure- 
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FIG. 1: (color online) (a) The single-apparatus model of re- 
peated measurements. All the measurements are performed 
with the same apparatus A. Before every measurement the 
state of A is initialized to a given state |app). (b) Multi- 
apparatus model for repeated measurements. Each measure- 
ment is done with a individual apparatus with initial state 
|app>. 



The dynamical model of a single quantum 
non-demolition measurement 



According to the dynamical theory of quantum non- 
demolition measurements pTjj . the measurement process 
can be described by the coupling between the to-be- 
measured system Q and the apparatus A. The total 
Hamiltonian of Q and A has an expression of conditional 
dynamics 



H M = J2\j)(j\®Hj 



(1) 



where Hj is the Hamiltonian of the apparatus A with 
respect to the j-th eigenstate | j) of the observable of the 
open system. Before the measurement, the system Q 
can be in any superposition state TZj—i Cj |j) , while the 
apparatus is set in a pure state |app). If the duration 
time of the measurement is tm, the measurement leads 
to the transformation 



^Q|j)|a PP )^^Q|. ? )|A 3 ) 



(2) 



where 



ment. In the multi-apparatus model, each measurement 
is achieved by an individual apparatus [36|. Namely, in 
every measurement, only one apparatus interacts with 
the to-be-measured system while all the others are left 
alone. In the recent experimental realization of QZE 
in cavity QED [Io[, in every measurement the state of 
the cavity field is measured by an individual ensemble of 
cold atoms. It can be considered as an illustration of the 
multi-apparatus model. 



In the current section, we prove that, regarding to the 
evolution of the to-be-measured system (in our case, it 
is the TLS together with the heat bath), the two models 
are equivalent. They can lead to the same evolution of 
the density matrix of the to-be-measured system. For 
the convenience of our calculation, in this paper we will 
use the multi-apparatus model in our discussion for QZE 
and QAZE from the next section. 



In the following we will first show the dynamical model 
of a single quantum non-demolition measurement, and 
then do the formal calculation for the evolution of the to- 
be-measured system under repeated measurements with 
single-apparatus and multi-apparatus model. Our result 
shows that the density matrix of the to-be-measured sys- 
tem has the same time evolution in both models. For 
the generality of our discussion, in this section the to-be- 
measured system is assumed to be a general multi-level 
quantum system. From next section we will focus on the 
system of TLS together with the heat bath. 



|Aj> = 



_ p -iHjTM 



app) 



(3) 



is the finial state of the apparatus with respect to the 
state \ j) of the system Q. Therefore, in the more general 
case, if the initial density matrix of the system Q is pq, 
after the measurement the density matrix of Q becomes 



(4) 



Then the effect of the quantum non-demolition measure- 
ment on the to-be-measured system Q can be described 
by the relevant decoherence factors (Ai\Aj) in the defi- 
nition (Fj} of the super-operator C. 



B. The single-apparatus model of repeated 
measurements 

Now we consider the case of repeated measurements. 
As shown in Fig. 1(a) and 1(b), we assume the measure- 
ment is performed once in every time region [^_i,£n] 
(n = 1, 2, t' = 0) with duration tm- During the mea- 
surement, the Hamiltonian of the system Q and the rel- 
evant apparatus is Hm- The length of the time intervals 
between every two neighbor measurements is assumed to 
be r. The system evolves freely with Hamiltonian Hq in 
the time between the measurements. 

In the single-apparatus model, all the measurements 
are performed via the same apparatus A (Fig. 1(a)). 
As we shown above, we assume the state of the appa- 
ratus A is initialized to a given state |app) before every 
measurement, and the information obtained by A from 
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the last measurement is "erased". This initialization or 
erasing process can be done via switching on the inter- 
action between A and an external reservoir R, like the 
spontaneous emission of the two-level atom. The similar 
technique has also been used in our bang-bang cooling 
scheme for the nano-mechanical resonator [46j]. 

Now we consider the evolution of the density matrix 
of Q in the repeated-measurement process. At the initial 
time t , the density matrix of Q and A is 



Pqa (t' ) = p Q |app)(app|. 



(5) 



During the time region from t' to ti, the first measure- 
ment is performed via the interaction between the system 
Q and the apparatus A while the interaction between A 
and the reservoir R is switched off. As shown above, at 
the ending time t\ of the first measurement, the density 
matrix of the system Q becomes 



PQ (t 1 )=£[ PQ (f )] 



(6) 



where the super-operator C is defined in Q. 

In the time region between t\ and t\, the Q-A inter- 
action is switched off. The system Q experiences a free 
evolution governed by the free Hamiltonian Hq, while 
the apparatus A and the reservoir R experiences a cou- 
pling which induces the initialization of the state of A. 
We denote the evolution operator given by the A-R cou- 
pling as uar- Then we have the total density matrix px 
of Q, A and R at time t[ as 



p T (t'l) = UpUARpT <4 



AR^F 



with 



up = cxp [—iHqt] 



(7) 



(8) 



It is pointed out that, since up and uar are the operators 
for different systems, we have 



[u Fi u A r] = 0. 



(9) 



Then we have the density matrix of the system Q at time 



PQ (t'l) 



Tv A r 

ufTtar 

ufPq (ti)u F 
UC [ PQ (t )] 



U F U A RPT (tl)u\ R u\ 



uarPt (h)u 



t 

AR 



(10) 



with the super-operator U defined as U [-••]= up [...] u F . 
On the other hand, the density matrix of the apparatus 
A at time t\ is given by 



Pa (t[) 



Tr 



QR 



uarUfPt {ti)u F u ! 



Tr, 



Tr fl 



UAR^Q 



F^AR 
Upp T (tl) U F 



AR 



UARPAR (h) U' AR 



|app)(app| 



(11) 



with the density matrix par (£i) of A and R at time t\ : 
p A R(h) = Tr Q [p T (t 1 )]. (12) 
In the last step, we have used the fact 



Tr, 



uarPar {h) u\ R = | app) (app 



(13) 



Namely, the A-R coupling can make the density matrix of 
A to decay to the unique steady state | app) (app |, which 
is independent on the density matrix of A at the time 
ti before the switching on of the A-R coupling. That 
"spontaneous-emission-like process" is the physical ex- 
planation of the "initialization" of the apparatus state. 
Eq. (fT3|) is applicable when the influence of the A-R cou- 
pling on the reservoir R is negligible in the total evolution 
from time t' to t' N . Then the density matrix of Q and 
A at the beginning time t[ of the second measurement 
becomes 



Pqa = Pq {t[) | app) (app |, 



(14) 



which has the same form as the one Pqa (to) at the begin- 
ning of the first measurement. Therefore we can straight- 
forwardly generalize our above discussion to the time af- 
ter t\. Finally we have the density matrix of the system 
Q at the time t' n after n measurements and n free evolu- 
tions: 



P Q (t' n ) = {UCT[pQ (t )] 



(15) 



C. The multi-apparatus model of repeated 
measurements 

Now we consider the multi-apparatus model for the re- 
peated measurements. In this model we assume there are 
many individual apparatus A^\ A^ 1 ', each of which is 
in the same state |app) before the measurements. In the 
m-th measurement in the time region [t^-i: tm\ j the sys- 
tem Q only interacts with the m-th apparatus A^ m \ and 
leave the other apparatus alone. In the time interval be- 
tween two measurements, the evolution of Q is also gov- 
erned by the same Hamiltonian Hq. Nevertheless, the 
initialization of the state of the apparatus is not required 
in this multi- apparatus model (Fig. 1(b)). 

At the beginning time t' of the first measurement, the 
density matrix of Q and A^ is 



Pqaw (t'o) = PQ (t'o) |app)(app| 



(16) 



At the time t\, after the first measurement, the density 
matrix of the system Q is also C [pq (t' )}. At time t[, the 
density matrix of Q is 



PQ (t[) = u F pQ (ti) up = UC [p Q (t' 



(17) 



Therefore, at the beginning time t[ of the second mea- 
surement, the density matrix of Q and the apparatus A^ 
is 



Pqaw (t'i) = Pq (ti) |app)(app| 



(18) 



with the same form as Pqaw (^o)- Then at time t' n after 
the n-th measurements we also have the density matrix 
of Q 

PQ (4) = {UC) n [p Q (f )] . (19) 

That is the same as the one (fT5|) given by the single- 
apparatus model. 

Therefore, the single-apparatus and multi-apparatus 
model for the repeated measurements can lead to the 
same evolution (fT9| of the to-be- measured system Q. If 
we only consider the evolution of the system Q under 
the repeated measurements, we can use either single- 
apparatus model or multi-apparatus model in our cal- 
culation. In our discussion in the following sections, the 
system Q includes the TLS and the heat bath which is 
coupled to the TLS in the time region of free evolution. 
We will use the multi-apparatus model to describe the re- 
peated measurements periodically performed on the TLS. 

III. REPEATED MEASUREMENTS ABOUT 
TWO LEVEL SYSTEM 

A. System and measurements 

In the above section we post our multi-apparatus 
model for the repeated quantum measurements. From 
this section, we show our dynamical approach for the 
QZE and QAZE of a dissipative TLS under repeated 
measurements. 

We consider a TLS coupled with a heat bath which is 
described multi-mode bosonic field. The Hamilto- 
nian of the total system is 

Hp = uj eg \e) {s) (e\+J2^k a l a k 
fc 

+ |e) (5) (ff| £ (g%a % + ,gl4) + h.c. (20) 

Here \g)^ and |e)( s ) are the ground and excited states 
of the TLS, aX and ag are respectively the creation and 

annihilation operators of the boson in the k-th heat-bath 
mode with frequency w^, while is the relevant coupling 
intensity between the boson and the TLS. In our Hamil- 
tonian Hp, the rotating wave approximation is not used 
so that the possible effects given by the counter-rotating 
terms can be included [13, 13 . 

In this paper we consider the evolution of such a dis- 
sipative TLS under repeated measurements. As shown 
in the above section, we use the multi-apparatus model 
to describe the repeated measurements. We also assume 
the n-th measurement is performed in the time region 
[*n-ij*n] ( n = 1j2>— ; t'o = 0) with the duration time 
tm (Fig. 2(a) and Fig. 2(b)). The time intervals between 
two neighbor measurements are also assumed to have the 
same length r. In these time intervals between the mea- 
surements, the system evolves freely under Hamiltonian 
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FIG. 2: (color online) (a) The multi-apparatus model of the 
repeated quantum measurements of a dissipative TLS. In each 
measurement, the TLS is coupled with an individual appara- 
tus and leaves other apparatus alone. The n-th measurement 
occurs in the time region between t' n _ 1 and t n . The TLS is 
also coupled with the heat bath, (b) The time sequence of 
repeated measurements and free evolution governed by Hf 
defined in (|20[) . (c) The definitions of the functions T n (t) and 
A n (t). 



Hp. In the multi-apparatus model, we assume there are 
many apparatuses A^\...,A^ n \... which can be individ- 
ually coupled with the TLS and distinguish between the 
states |e)( s ' and \g)^ . In this paper we denote \)^ 
for the quantum state of the TLS, |)( B ) for the state of 
the heat bath, |)( n ) for the state of the n-th apparatus 
and |)^ for the state of all the apparatus. Before 
the measurement, every apparatus is initially in a pure 
state. During the n-th measurement, the TLS is coupled 
with the n-th apparatus A^ and decoupled with all the 
other apparatuses ^4( m ^™). 

As shown in Sec. II, in the quantum non-demolition 
measurements the coupling between the system and the 
apparatus A^ has an expression of conditional dynam- 
ics: 

= |e) (S) (e| ® ffW + | 5 ) (S) (.9l ® H<f\ (21) 
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where (H^) is the Hamiltonian of the apparatus 

A^> with respect to the state |e)( s ) (\g) {s) ) of the TLS. 
In this paper, we assume the duration time tm of the 
measurement is so small that the interaction between the 
TLS and the heat bath and the decay of the excited state 
can be neglected during the measurement. 

In the n-th measurement, if the state of the TLS before 
the measurement is a|e)( s ) + /3\g)( s \ and the one of the 
apparatus A n is |app)(™', then the transformation given 
by the measurement can be described as 

(a|ep+/%)^) |app)M 

'™|e)< s )|A e ) (n) + P\g} (S) \A g ) {n) 



— > ae 



with \A g )^ and \A e )^ the states of the apparatus at- 
tached to the ground and excited states of the TLS 



\A g ,e) in) 



|app) 



(») 



(22) 



The ideality of the measurement is described by the over- 
lap of the states \A g )^ and |A e ) < -"' 1 or the decoherence 
factor 



7ne 



(n) 



{A e \A g ) {n) ] 7„,6» n e Reals. 



(23) 



In the case of ideal projective measurement, these two 
states are orthogonal to each other and we have -f n = 0. 



Here |^ (t)) is the state of the total system in the 
Schrodingcr picture and Hq is given by 

H (t) = u; eg \e)^(e\ + J^w^ + Y, H M r n (*) (28) 



In the interaction picture, the quantum state |\& (t))j sat- 
isfies the Schrddinger equation 



i-\*{t)) I = H I (t)\9(t)) I 



(29) 



with the Hamiltonian Hj (t) given by 

H I (t) = \e)W(g\f B (t)f A (t) + h.c, (30) 
where the operators fs (t) and Ja (t) are defined as 

hit) = E(^ e_iA£ * + ^4 e<Wes+ws> )' (31) 



f A (t) = J2 A n(t)M n . 



n=l 



Here the detuning Ar takes the form 



(32) 



(33) 



and the unitary operator M n is given by 



B. The total Hamiltonian in the interaction picture 

Under our above assumptions, the quantum evolution 
of the total system including the TLS, the heat bath 
and the apparatuses is governed by the time-dependent 
Hamiltonian 

H (t) = u eg \e)^ (e\ + Wjs4°iS + E H m T " (*) 



+ (\e)^(g\ + \g)^(e\) i^g^ + h.c.j ^A n (t), 

(24) 

where the functions T n (t) and A„ (t) are defined as in 
Fig. 2(c): 



and 



r„ (t) 



A„ (t) 



0; otherwise 



Me[fn,fJ 

0; otherwise 



(25) 



(26) 



To solve the quantum evolution of our system, we use 
the interaction picture where the quantum state |\t (t))i 
is defined as 



|* (i))/ = exp 



i / H (t')dt' 
Jo 



!*(*)>• 



(27) 



M n = Y[ exp 



exp 



(34) 



i=i 



With the help of the interaction picture, the effect of 
the measurements is packaged in the definition of the 
operator M n . As shown below, in this interaction picture 
our calculations are significantly simplified and we can 
express all the effects from the measurements in terms of 
the decoherence factor 7„e* e " defined in Eq. (|23| . In the 
following two sections we will derive the short-time and 
long-time evolutions of the TLS with the calculations in 
this interaction picture. 



IV. SHORT-TIME EVOLUTION: FIRST-ORDER 
PERTURBATION THEORY 

In this section we calculate the short-time decay rate 
of the excited state |e)( s ) of the TLS under repeated 
measurements. For simplicity, we only consider the zero- 
temperature case where the initial state of our system at 
t = t' = as 



|*(0)) / = |e)( s )|vac)( s )ni a PP) (m) ' 



(35) 



where |vac)( s ) is the vacuum state of the bosonic field 
and |app)( m ) the initial state of the m-th apparatus. We 
consider the evolution of the system from t = to a finial 
time t = tp. which, for simplicity is assumed to be an 
integer multiple of the period t+tm of the measurements. 
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The finial state (£f)) / of the total system is the solu- 
tion of the Schrodingcr equation (|29j) . It can be expressed 
as a functional series of Hj(t): 

I* (tf )>/ = U~i JJ H I (t)dt + |* (0)} 7 . (36) 

The survival probability P e (t F ) of the state |e)^ s ^ is given 
by 

P e {t F ) = 1 i(*(t P )\(\g)M(g\)\*fo)h- (37) 

We can define the short-time decay rate R(t, t F ) of |e)' s ' 
as 



P e (t F ) = l-R(T,t F )t F . 



(38) 



Therefore, R(r,t F ) can also be expressed as a functional 
series of Hi (t) . When the total evolution time t F is small 
enough, we can only keep the lowest-order term (in our 
problem it is the second-order term) of Hj(t) in R(r,t F ). 
Then we have the short-time decay rate 



R{tM^ T Tf Y1 \9k\ 2 ^ 2 (Agr/2) 

k 



N 



N TO— 1 

2 EE Rc 

m— 1 n—1 



^i(n-7n)A^(r+T M ) 



n 



l=n+l 



(39) 



where 



iV 



*f 



tm 



(40) 



In the case of ideal projective measurements with neg- 
ligible duration time or tm = 0, we have 7; = and then 
the short-time decay rate R(r, t F ) in Eq. (|3"9"|) becomes 
a ^-independent one: 

R(t, t F ) = R plo (t) = r £ lysine 2 (Agr/2) , (41) 



It means that, in the cases of non-ideal measurements 
with 7 max < 1, the QZE also occurs in the limit that the 
measurement are frequent enough. 

In Refs. [34l - [37| . the authors have calculated the total 
survival probability P e (t F ) as the product of the one 
after each measurement, i.e., the relationship 



P(t F ) = P{t) 



N 



(44) 



is assumed. This simplification yields that the short- 
time decay rate R(t, t F ) is the one in a single period of 
measurement, or 



R(r,t F ) =R(t,t). 



(45) 



In our case of t ^ tm, Eq. (|45|) leads to R(r,tp) = 
R-pro (t). Therefore, the intuitive treatments (|44[) and 
(|45p arc reasonable in the cases of ideal measurements 
with 7„ = 0, and required to be improved to Eq. (|39p for 
the cases of non-ideal measurements and non-zero 7„. 



A. The effect given by non-ideal measurements 

To further explore the physical meaning of the short- 
time decay rate given by Eq. (|39[) . wc consider the simple 
case with identical non-ideal measurements. In this case 
the decoherence factor 7„e* e " takes n-independent value 
je ie (7 < 1). In the large N limit with TV > 1/(1 - 7), 
Eq. ([3"9"1> gives a ^-independent short-time decay rate: 



R(r,t F ) w i? mca (T) 



dt]G (5 + u eg ) h( 7 ,e- 77) sinc 2 |.(46) 



Here we have used the approximation t/(t + tm) ~ 1- In 
the above expression the spectrum function G (w) of the 
heat bath defined as 



(47) 



which is the result given by Ref. with the postula- 
tion of wave function collapse. Therefore, the QZE or 
QAZE based on the ideal decay rate i? pro (t) can also be 
obtained with our model. 

In the case of periodic non-ideal measurements, if the 
common upper limit of the modului j n of the decoherence 
factors in the measurements are smaller than unit, i.e., 
we 7„ < 7max < 1 for any n, then Eq. ([39j) gives 



R(r,t F )<rJ2\9k\ 

k 



2jn 



l-7n 



nc 2 f — ^ ) (42) 



Therefore, for a fixed value of the total evolution time 
t F , we have 



lim R(t, tp) = 0. 

T— >0 



(43) 



and the function h (7, x) is given by 

1-7 2 



h (7, a) 



1 + 7 2 — 27 cos x 



(48) 



It is obvious that, in the case 7 = we have h (0, x) = 1 
and the decay rate i? m ca (t) in Eq. (j4"6")l returns to the re- 
sult i? pro (r) given by ideal projective measurements. In 
the case of repeated non-ideal measurements, the func- 
tion h (7, 9 — 7]) in the right hand side of Eq. (|46|) would 
tune the shape of the function to be integrated and then 
change the value of i? mca (t) . 

In Fig. 3 we plot the function h (7, 6 — rj) with respect 
to different values of decoherence factor je 1 ® . It is easy 
to prove that h (7, 6 — rj) takes peak value at the points 
i] = 6 + 2mr (n = 0, ±1,....) with the width of the order 
of arccos7. Therefore, the effects given by h (7, 9 — rj) on 
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FIG. 3: (color online) The function /i(7, 8 — rf) defined in 
Eq. (|48p . Here we plot the curves with 7 = 0.3, 8 — (green 
line with open circle); 7 = 0.5, 8 — (pink dashed line); 
7 = 0.8, 8 = (black solid line) and 7 = 0.8, 8 = 2tt/3 (red 
dotted line) . It is clearly shown that the position of the peak 
of /i(7, 8 — tj) is determined by the value of 9, while the width 
of the peak decreases when 7 is increased. As a comparison, 
we also plot the function sine 2 (77/2) (blue dashed dotted line). 



the decay rate i? mea (r) in Eq. (|46|) seriously depend on 
the values of both 7 and 8. 

To illustrate these effects, we calculate the decay rate 
■Rmea (t) of a TLS in an environment with noise spectrum 

G(w) = G H w) = ^ 0+(w) ) , (49) 

[ 0; cj<0 

which is, up to a global factor, the same as the noise 
spectrum of the 2p-ls transition of the hydrogen atom 
[49l l5(ij . As in the case of realistic hydrogen atom, here 
we also take ui c = 549. 5w e9 . The decay rates i? mca (t) 
with = 0, 7r/2, 7r, 37r/2 and various values of 7 are shown 
in Fig. 4. 

Our results show that, with any values of 7 and 6, the 
short-time decay rate i? m ca (t) approaches to zero in the 
limit r —> 0. That is, the QZE always occurs when the 
measurements are repeated frequently enough. This is 
consistent with the conclusion in Eq. (|42|) . On the other 
hand, in the limit r — > 00, i? mea {t) approaches to the 
natural decay rate Rgr given by the Fermi golden rule: 



R gr = 2ttG (w efl ). 



(50) 



Therefore, the quantum Zeno and anti-Zeno effects only 
occur when r is small enough. 

When the phase 8 = 0, as shown in Fig. 4(a), both 
the QZE and QAZE can also appear even when the de- 
coherence factor 7 is nonzero. However, when the value 
of 7 or non-ideality of the measurements is increased, 
the region of r for the appearance of QZE and QAZE 
becomes narrower. That is, in the case of large r, the 
QZE and QAZE with large 7 are always less significant 
than the ones with small 7. In this sense the two effects 
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FIG. 4: (color online) The short-time decay rate -R me a(T) de- 
fined in Eq. (|46p of a dissipative TLS under repeated mea- 
surements. Here we take the noise spectrum in Eq. (|49[) and 
illustrate the cases with 8 = (a), 6 = 7r/2 (b), 9 = 7r (c), 
# = 3tt/2 (d) and 7 = (blue solid line), 7 = 0.3 (green 
dashed line), 7 = 0.8 (black dashed-dotted line). We also 
plot i? mca (r) (red open circle) defined in Eq. ( I5II) given by 
frequent phase modulation pulses with approximation (|53l) . 



are suppressed by the non-ideal measurements with real 
decoherence factors. 

For the cases of non-zero phase 8, the effects given 
by the non-ideal measurements are from both the non- 
zero value of 7 and the repeated phase modulation given 
by the phase factor e . As shown in Figs. 4(b)-4(d), 
the final behavior of i? mC a sensitively dependents on 8. 
Especially, with specific values of 8 and 7, the non-ideal 
measurements can enhance either QZE or QAZE. Both 
the r-region where the QZE occurs and the peak value 
of the decay rate in the QAZE can be enlarged by the 
complex value of the decoherence factor -fe 10 . Finally, the 
total region of r for the occurrence of QZE and QAZE 
with non-zero 7 is possible to be the same as the one 
for the cases of 7 = 0, or even larger than the latter. 
That is due to the complicated behavior of the function 
h(j,8-r)). 



B. The quantum Zeno and anti-Zeno effect given 
by phase modulation pulses 

As shown in the above subsection, in the cases of non- 
ideal measurements with 7 7^ 0, both the QZE and QAZE 
are possible to be enhanced in the cases with nonzero 
phase shift 8 in each measurement. Actually in the most 
extreme cases of 7 = 1, the two effects can also appear 
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with nonzero phase 9. In that case, the measurements 
reduce to the phase modulation pulses which can induce 
periodical jumps for the relative phase between \e)^ s > and 
\g)( s) . The QZE and QAZE given by the phase modula- 
tion pulses have been discussed in Ref. [Io, in detail. 
Here we show that our dynamical model is also able to 
describe these effects. 

If the pulse is repeatedly performed with period r and 
negligible duration time tm, then the short-time decay 
rate in Eq. ([59"]) becomes 

R(r,t F ) = i? PM p(r) 

= J d V G (J + uj eg ) g {9 - r,) sine 2 (|) ,(51) 



where the function g(x) is defined as 
1 sin 2 (Nx/2) 



9 0) 



N sm 2 (x/2) 



(52) 



This result is essentially equivalent to the one in Refs. [40|, 

m. 

From now on we assume the function G{uj + uj eg ) has a 
finite width Aw, i.e., G(uj + uj eg ) takes nonzero value only 
in the region —Auj < u> < Auj. Under this assumption, 
when the evolution time tp is large enough so that 2-k /tp 
is much smaller than the width Aw, we have 



g(x) « 2k Hx + 2im) . 



(53) 



If the period r of the pulses is short enough so that r < 
k j Auj and 



G(-7r/T + UJeg) = G(k/t + UJ eg ) ~ 0, 



(54) 



the function G(t]/t + uj eg ) is localized in the region — ir < 
•q < 7r. Then the simplification (|53[) implies 



-Rpmp ~ 2nG \ - + uj eg I smc ^ 



-2-kG 



9 2-k 



+ U! eg ) smc 



7T ,(55) 



where we have assumed 9 € [0,2-7r]. In this 
result of Eq. (|54| , one can further find some special 9 so 
that 9/t — 2tt/t < —Aus and 9/t > Auj, which makes 
G(6/t + UJeg) - G{9/t - 2tt /r + uj eg ) - or i? PM p ~ 
0. There are also other special angles 9* which makes 
G(9* /r+ujeg) or G(9*/t — 2-k / r+uieg) take the maximum 
value G max of G(uj). Therefore, the decay rate i?pMP 
can be tuned in the broad region between and some 
maximum value which is of the order of 27rG max . This 
tuning effect is also predicted in Refs. fifl l4l|. 

We also point out that, in the limit 7 — >■ 1 the function 
h (7, x) in the last subsection has the same behavior 



h (7 -t l,x) w 2tt 



00 

E 



2n?r) 



(56) 



as g(x) in the large N limit. Therefore the results in 
Eq. ([51) for 2n/tp << Auj can also be obtained from 
Eq. (0SD. 

In Fig. 4 we also plot the short-time decay rate i?pMP 
under approximation (|53[) with the spectrum (|49|) and 
different values of the angle 9. It is shown that, when 7 
is close to unit the short-time decay rate i? mea in Eq. (|46|) 
is quite close to -Rpmp ■ In this case the behavior of the 
short-time decay rate is dominated by the repeated phase 
modulation. 



V. THE LONG-TIME EVOLUTION: RATE 
EQUATION 

In the above section we have considered the short-time 
evolution of the dissipative TLS under repeated mea- 
surements or phase modulation pulses. We obtain the 
short-time decay rates via perturbative calculation based 
on our pure dynamical model of repeated measurements. 
The perturbative approach is simple and straightforward. 
Nevertheless, the results arc only applicable when the to- 
tal evolution time tp is short. 

In this section, we go beyond the short-time calculation 
and consider the long-time evolution of the TLS under 
repeated measurements. The problem wit h proje ctive 
measurements has been considered in Refs. |42h44| in a 
semi-classical approach with measurements described by 
a stochastic term in the Hamiltonian. Here we provide 
a full-quantum theory which can be used for the cases 
of either ideal or non-ideal measurements. The previous 
results [10, SU on the long-time evolution of a dissipative 
TLS under repeated phase modulation pulses can also be 
derived in our approach. 

For simplicity, here we assume the measurements are 
identical with the decoherence factor je 10 . We first de- 
duce the general form of the rate equation of the TLS 
in terms of the effective time-correlation function of the 
environment, and then derive the simplified form of the 
rate equation under the time-local and coarse-grained ap- 
proximation. 



A. The general rate equation and effective 
time-correlation function of the environment 

To derive the rate equation for the TLS, we first con- 
sider the TLS and the apparatuses as a total system in- 
teracting with the heat bath. The evolution of the den- 
sity matrix p( SA ^ of such a combined system can be de- 
scribed by master equation given by Born approximation 
(Eq. (9.26) of Ref. [EU): 



= — I cIsTib 
Jo 



H z (t), iJ 7 (s),P (SA) (s)®P (B) (57) 
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where Hi (t) is defined in Eq. (|30|) and is the initial 
density matrix of the heat bath, which is assumed to be 
in the thermal equilibrium state at temperature T. It 
is pointed out that, in Eq. (|57[) we do not perform the 
Markovian approximation. 

The evolution equation for the density matrix p^ (t) 
of the TLS can be obtained by tracing out the states of 
the apparatuses in Eq. (f57|) . According to Eq. (|30|) , in 
the calculation we need to evaluate the values of 



Tr A 



f A {t)P A { S ) P {SA) (a) 
f A {s)f A {t)p {SA) («) 
! A [t)fA{ S ) P {SA) (a) 
f\(s)f A (t)p iSA) (a) 



(58) 
(59) 
(60) 
(61) 



with s <t. Noting that the product of f A (t) (/t (t)) and 
f A {s) (f A {s)) is nonzero only when there exist integers 
rit and n s so that A„ t (t) = A„ s (s) = 1 or 

3 n t ,n s , so that t E [t' nt ,t nt ] ; s 6 [4 S >*«»] ■ (62) 

When this condition is satisfied, we separate all the ap- 
paratuses into two groups: 

X : {the apparatuses A\, A 2 ..., A Ub } 
Y : {the apparatuses A ns +\, A 7ls+ 2, ...} 

Obviously, X is the group of the apparatuses which inter- 
act with the system before the time s, while Y includes 
the ones interact with the system after the time s. There- 
fore the density matrix p( 3A > (s) can be written as 



(63) 



where p( sx *> (s) is the reduced density matrix of the TLS 
together with the apparatuses in the group X, and p\ is 
given by 



|a P p)( m >(app|. 



(64) 



m—n s -\-l 



On the other hand, Eqs. ([3l|l. (|32jl and ([3"3]) „ together 
with the condition s < t, yield 



fA(t)f A (a) 



n ex p 

i=n,+l 



MP 



exp 



(65) 



Namely, the operator f A (t) f' A (s) only operates on the 
state of the apparatus in the group Y. Then the quantity 
in Eq. (|58|) can be obtained as 



Tr A [f A (t)f A (s)p( SA Hs) 
p^ (s) Tt a \f A (t)f A (a)p* 



(66) 



It can be calculated easily with the simple form in 
Eq. ([64]l of Pq . The terms in Eqs. (|59ti61|) can be evalu- 
ated in the similar approach. 
Then we get the rate equation 



I™ 



daF A + B ) (t,a)P e (a) 
dsF^(t,a)P g (a) 



(67) 



where P e (i) =< s ) (e^l {t)\e)^ and P g (t) = 1 - P e (t) 
are the probabilities of the excited and ground states 
of the TLS. In the cases of Sec. IV where the initial 
state of the TLS is assumed to be |e), the P e (t) de- 
fined here becomes the survival probability defined in 
([37]l . In Eq. ([67]l . the effective time-correlation functions 

Fab (*> s ) °f ^ ne environment are defined as 



Fab (*> s ) = 2Re 9b' (* - s ) 9a (t, s) 



(±) 



(68) 



where the bare correlation- functions of the heat bath 
are given by 

9 i B +) (t-s)= r Tr B [f B (t)fUs)p {B) ] (69) 
= J2 \9k\ 2 \( n k + 1) e~ tA ^- s) + n% e<^ +u "o)^ 



g ( B - ) {t-s)=Tt B P B (t)f B {s)p 



.(B) 



(70) 



(*-«) -|- + l) e -' ; ("£+"eg)(t-s) 



with rig = Ti's[ata- the average number of the bo- 
son in the k-t\i mode of the heat bath. The correlation- 
function of the measurements is defined as 



g A (t, a) = Tr A f A (t) f\ (s) pi 



(71) 



It is easy to prove that, the functions g A (t, s) and 
g^ (t — s) decrease when the absolute value of \t — s\ 
increases. When t and s satisfy the condition f|62[) we 
have 



g A ^ s ) =7 ^ t -n s ) e i(n t -n s )6 



(72) 



If the condition ([62]) is violated we have g A (i, s) = 0. 

The rate equation (|67|) shows that, the evolution of the 
probabilities P 6j9 (s) of the excited and ground states of 
the TLS is governed by the time-correlation functions 
F AB (*j s )j which are given by both the time-correlation 
function of the heat bath and the decoherence factors of 
the measurements. The measurements tune the corre- 
lation function F AB (t, s) through the function g A (t, s). 

Especially, the trail of F A J s ) m ^ ne long-time-interval 
region with large \t — s\ would be suppressed by the fac- 
tor <-y("*- n =) in the function g A (t, s) defined in Eq. (J72J). 
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FIG. 5: (color online) The effective time-correlation function 

Fab (t, 0) defined in Eq. (|68[l of the environment of a dissipa- 
tive TLS under repeated measurements. Here the heat bath is 
assumed to be at zero temperature with Ohmic spectrum (]73[l , 



where uj c = 500oj eg . The incoherence factor of the non- ideal 
measurements is 7 = 0.5 and 9 = 0. We plot the behaviors of 
with the measurement period r = 0.003/cj es (blue line 
with open diamond), 0.005/u>e 9 (red line with open square), 
0.01 /uj eg (black line with open triangle) and 0.1/u) eg (green 

line). It is clearly shown that the long-time trail of F A + g{t, 0) 
is suppressed by the frequent measurements with small r. We 
also take the limit tm = 0. 



To illustrate the effects given by the repeated measure- 
ments to the effective correlation function F^g (t,s), in 

Fig. 5 we plot Fj^g (t, 0) for a TLS in a zero-temperature 
environment with Ohmic spectrum 



G(u) = G (u) = we-"/" 



(73) 



It is clearly shown that when the increasing of the fre- 
quency of the measurements, or the decreasing of the 
time interval r between measurements, can lead to the 
suppression of the long-time trail of F A ^g (t, 0). 

One important parameter for the correlation function 
Fj^g (t, s) is the effective correlation time tf, which gives 

F { ab (t, s) ~ for \t-s\> TF- If tf is small enough so 
that the variation of the probabilities P e (t) and P g (t) 
is negligible in the time interval Tp, we can perform the 
time-local approximation 

P e , g (s) « P e , g (t) (74) 

in Eq. (|67p and significantly simplify the rate equation. 
Eq. (Ell) yields 

t f = min {ta,t b } , (75) 

where ta and tb are the correlation times of the functions 
qa (x) and ■ From Eq. (|T2j) we have 

t a ~ (1 - 1/ In 7) (t + t u ). (76) 

In the case of ideal projective measurements we have 
Ta = t T Ty[ . Therefore, Eq. ([75]) shows that frequent 
measurements with small period t can reduce the effec- 
tive correlation-time of the environment experienced by 
the dissipative TLS. 




t (units of 1/ro ) 
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FIG. 6: (color online) The decay rate R e (t) in Eq. (JTHJ. Here 
we take the measurement period r = 0.1/cj eg . Other parame- 
ters are the same as the ones in Fig. 5. As proved in Eq. (|81|) . 
in the long-time limit the decay rate is a periodic function of 
t. (Here we also take the limit tm = 0.) The discontinuity 
of the function R e (t) at t = nr (n = 1, 2, 3, ...) is due to the 
jumping behavior of the function gA (t, s) in Eq. 1)721) . 



B. The master equation under time-local and 
coarse-grained approximation 

In the following we assume the effective correlation 
time Tp is small enough and the time-local approxima- 
tion (|T4"|) can be used. Then the rate equation (|6"T|) can 
be simplified as 



dt 



P e (t) = ~R e (i) P e (t) + Rg (t) Pg (t) , 



(77) 



where the time-dependent decay rates R e (t) and R g (t) 
are given by 



Re,g (t) 



(78) 



On the other hand, Eq. (|72|) implies that 

g A (t, s) = g A {t + t + t m ,s + t + tm) ■ (79) 
This relationship, together with the definitions of 



F 



AB (M): S iveS 



Fab (*> s ) = f ab ] (t + r + t m ,s + t + t m ) 



(80) 



Therefore, when t is much longer than the effective cor- 
relation time Tp, we have 



Re,g (t) 



f : dsF$(t,s) 
Jo 
ft 

/ dsF^B {t + T + T M ,S + T + T M ) 

Jo 



t+T+TM 



dsF^ {t + T + TM, S) 



T+TM 

t+T+TM 



/■t+T+TM 

/ dsFX% (t + T + TM, S) 

Jo 



Re,g (t + T + TM), 



(81) 
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and then the decay rates become periodic functions of t 
with period r + tm , which is the same as the period of 
the measurements. In Fig. 6 we plot the R e (t) for the 
system in the calculation of Fig. 5. The periodic behavior 
of R eg (t) in the large t case is illustrated clearly. 

If the measurements are frequent enough so that the 
variation of the probabilities P e _ g (t) in the time interval 
t + tm can be neglected, we can further perform the 
coarse-grained approximation and obtain the Markovian 
rate equation 



0.5 



P e (t) = -R GG P e (t) 



dt 



(t) 



(82) 



with the coarse-grained decay rates 



R 



CG 
e,9 



lim 

N^oc 



N(t + t m ) Jo 



N(t+t m ) 



R e , g (t) dt. (83) 



It is easy to prove that in the zero-temperature case we 
have 



R 



CG 



1 f f 

lim — / dt dsg^\t — s)gA(t, s) + h.c 
t.F J Q J 

9 ( B\t - s)g A (t,s) 



lim — / dt ds 

t F ^oo t F J J 



lim — 

tp^oo tp 



dfcfT z (t)|tt(0))j 



(84) 
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FIG. 7: (color online) The probability Pf of the excited 
state of the TLS in the steady-state. Here we plot P| given 
by Eq. I|88p. We take the noise spectrum in Eq. and 
illustrate the cases with (a) 9 = 0, (b) 9 = 7r/2, (c) 9 = n 
and (d) 9 = 37r/2, and 7 = (blue solid line), 7 = 0.3 (green 
dashed line), 7 = 0.8 (black dashed-dotted line). In (b-d) 
we also plot the probability P% (red open circle) given by 
frequent phase modulation pulses. 



with Hi{t) and |#(0))j defined in Eqs. fj30]) and ([35]). 
Therefore R^ G at zero temperature is reduced to the 
short-time decay rate R moa (T) being defined in Eq. (|46| . 
Namely, the decay rate given by calculation for the short- 
time evolution also governs the long-time evolution when 
the effective correlation time tf and the period of the 
measurements t + tm is short enough. 

Finally we discuss the behavior of the coarse-grained 
decay rate R GG of the ground state. It is well-known 
that, for a dissipative TLS without measurements, the 
decay rate of the ground state is usually negligible in 
the zero-temperature case. However, due to the counter- 
rotating terms of the Hamiltonian (f2TJl) , the correlation- 
functions defined in (|71[) and the ground-state decay 
rate are not absolutely zero. When the periodical mea- 
surements are performed, the value of the ground-state 
decay rate R CG is also varied by the measurements, and 
can take significant value even in the zero-temperature. 

For simplicity, we consider the simple cases with peri- 
odic identical measurements which has the decoherence 
factor rye 10 (7 < 1), i.e., the cases discussed in Sec. IV 
A. The straightforward calculation shows that in such a 
case the coarse-grained decay rates are given by 



fle,9M= / dr,G(^±uj eg ) h(rf,e- V ) sinc 2 |, (85) 

where "+" stands for R GG and "-" for R GG . In the 
above expression, r is length of the time interval between 



two measurements. The function h and the spectrum G 
of the heat bath are defined in Eq. (@8]) and (|47| . 

Therefore, when there is no measurements, or in the 
limit t — > co, we have the decay rates of ground state 
and excited state 

R gg (t -»• 00) « G (±u> eg ) J drjh (7, 6 - r,) sinc 2 |. (86) 

Since all the frequencies of the heat-bath are positive, we 
have G (— u) eg ) = and then R cg (t -)• 00) m 0. 

In the presence of periodic measurements, the decay 
rates R CG are given by the overlap of the spectrum 

G (i]/r ± Lo eg ) and the function h (7, 9 — rj) sine 2 (77/2), or, 
roughly speaking, given by the values of G (tj/t ±oj eg ) 
in the region r\ G [— it, tt]. Wc assume the function 
G{x) takes nonzero value in the region x G [0,f2], then 
G (t]/t ± to eg ) is nonzero only when 77 G [=Fw es r, T^egT + 
Qt]. Therefore, when the measurements are more fre- 
quent, or the time interval r becomes smaller, the non- 
zero region of G (tj/t — Lo eg ) has more and more overlaps 
with the region [— 7T, tt]. Then the decay rates R GG can 
be significant. In the limit r — > 0, both of the two func- 
tions G {t)/t ± (jj eg ) take non-zero values only in a small 
region around rj = 0. Then we have 

R GG (r -»• 0) = R cg (t -+())« r/»( 7l 6) [ d£G (£) . (87) 



The influence of the finite value of R GG can be observed 
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from the steady-state solution of the coarse-grained rate 
equation ([82jl . which describes the result of the long-time 
evolution of the dissipative TLS under repeated measure- 
ments. According to Eq. (|8"2"j) , in the steady state, popu- 
lation probabilities Pf and Pf of the TLS in the ground 
and excited states can be expressed as a function of the 
time interval r of the measurement: 

9 1 ) R GG {r) + R gg (t) ' e R GG (r) + R gg (t) ■ 

(88) 

Therefore, the non-zero decay rate R GG of the ground 
state leads to the non-zero population probability Pf' of 
the excited state. 

To illustrate the effects given by finite R^ G , in Fig. 
7(a)-7(d) we plot the population probabilities Pf with 
respect to different values of r and the decoherence factor 
je l6 . As in Fig. 4, we take the noise spectrum in Eq. (|49| . 
It is clearly shown that when t becomes small, the prob- 
ability Pf becomes non-zero. In the limit r — s> 0, Pf 
approaches 1/2 which implies R GG = R GG . In Fig. 7 we 
also plot Pf for the cases with phase modulation pulses 
rather than repeated measurements. The behavior of Pf 
is quit similar with the one given by repeated measure- 
ments. The non-zero decay rate R GG for the cases with 
phase modulation pulses is also obtained in Refs [13, EH . 



pulses can also be derived in our framework as a special 
case of non-ideal measurement. 

Based on our model, we derive the short-time decay 
rate (|39[) which implies that the QZE and QAZE may be 
enhanced by a non-ideal measurement with complex de- 
coherence factor. The long-time rate equation (j6"T|) is also 
obtained in terms of the effective time-correlation func- 
tion F^b , which describes the adjustment of the noise 
spectrum from the repeated measurements. The rate 
equation also shows that, the decay rate of the ground 
state of the TLS may be changed to non-zero value by 
the repeated measurements, and then the steady-state 
probabilities of the ground state and the excited state 
are also varied. 

The effects of non-ideal measurements on QZE and 
QAZE obtained in our model can be observed in the ex- 
periments where the system-apparatus interactions are 
well controlled. Such systems are possibly to be realized 
by nuclear magnetic resonance or solid-state quantum de- 
vices. Although we illustrate our model with a TLS, all 
the techniques in this paper, including the choice of the 
interaction picture in Sec. II. C, the time-dependent per- 
turbation theory in Sec. Ill, the master equation ([57| 
and the separation of the apparatuses in Sec. IV. A, can 
be used in the cases other than TLS. Therefore our model 
presented in this paper can also be straightforwardly gen- 
eralized to the discussions of QZE and QAZE of multi- 
level quantum system. 



VI. CONCLUSION AND DISCUSSION 

In this paper we provide a complete dynamical model 
for the evolution of dissipative TLS under repeated quan- 
tum non-demolition measurements. Our model gives an 
explanation of QZE and QAZE without the wave func- 
tion collapse postulation. The effects given by non-ideal 
measurements are naturally obtained in our model. The 
QZE and QAZE given by repeated phase modulation 
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